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I. INTRODUCTION 



The use of resonant (Gamow) states in nuclear structure calculations was proposed years ago by the Stockholm- 
Debrecen group 1^] and since then the notion has been widely applied to a variety of physical situations, with a 
remarkable success. The reader is kindly referred to Refs. [Q| j2| and references therein for a comprehensive introduction 
on the subject. Although Gamow's idea of decaying states was immediately recognized as a major breakthrough in 
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its use in modern nuclear structure calculations was delayed for nearly forty years until the 
shows that single-particle basis (Beggren's basis) can accommodate single particle resonant 
states of complex energy. Berggren's suggestions were adopted by Liotta and co-workers |^ thus given structural 
identity to a modern view of the continuum and its effects upon nuclear structure observables |2|. Mathematically 
speaking the Stockholm approach is, as a matter of fact, based on the identification of single-particle resonances, in 
the standard one-body nuclear central potential |^ , and in the numerical calculation of the needed matrix elements of 
nuclear two-body interactions [||. So far, the numerical treatment of resonant states, either by performing explicitly 
the needed integrals in the complex plane [Q Q or by projecting them on the real axis as non-overlapping states 
with a Breit-Wigner broadening, has proved its feasibility. A parallel formal development on Gamow's resonances 
followed from the work of Bohm and Gadella, which is documented in a rigorous mathematical way in |^]. However, 
a link between these two fronts of research on Gamow resonant states is missed. Particularly, we have noticed that 
basic notions, as the one of the value of the expectation value of the Hamiltonian on a resonant state, still need to be 
clarify and/or discuss in detail. In this paper we focus on the question about the definition of the expectation value 
of the Hamiltonian on a resonant state. We shall illustrate this point as follows. In ordinary Quantum Mechanics 
the mean value of the energy on a state represented by the density p is given by Tr pH, where H is the Hamiltonian. 
If p represents a pure state it has the form p = ^^'^ the quantity TrpH gives the usual formula for the mean 

value of the energy on the state which is written as {ip\H\ip). Gamow resonant states (or Gamow states, for 
short), are state vectors representing exponentially decaying states. They are not ordinary quantum states. Such a 
quantum state has a Breit-Wigner energy distribution, which is strictly nonzero for all values of the energy E, i.e. 
—oo < E < oo. This contradicts the fact that the spectrum of the physical Hamiltonian H should be lower-bound. 
Therefore, the Gamow state cannot be represented as a vector in a Hilbert space in which _ff is a self-adjoint operator 
§ (lO| |l| In addition, if (fi is a Gamow state, it must fulfill a condition of the kind e '^*^(p = Ae ''^ip, A being 
a phase {A = e'^^^^ with (j)(t) real). Thus, Hip — {<l){t)/t + ij)(p and the Gamow vectors must be eigenvectors of the 
total Hamiltonian with complex eigenvalues. Obviously this condition cannot be satisfy in the ordinary Hilbert space, 
since H must be self adjoint. In fact, Gamow states can be properly defined as functionals on a certain space of test 
vectors, in the same form that the generalized eigenstates of the position and momentum operators are defined [Tl| ] 
p2| . The space of functionals contains also the Hilbert space of ordinary states, so that no information is lost from 
ordinary Quantum Mechanics. Thus, we have a triplet of spaces ^ C H C , called the rigged Hilbert space or 
Gelfand triplet. 

Here Ti. is the Hilbert space, <I> the space of test vectors and "I>^ its antidualQ. In this context, Gamow states will 
exist in the dual space @] At this point a difficulty arises, namely: (E"^ is not an inner product space. If we 



£ if it is a mapping from $ into the complex plane C with the following conditions (i) Antilinearity: F{aip + I34>) = 
a*F{if) + /3*F{(j)), for all a, [3 £ C and all (/> £ $. (ii) Continuity: F must be a continuous mapping. is a topological 
vector space with a topology which is coarser than the Hilbert space topology on Ti. 
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represent the Gamow vector as |/o), the bracket (/o|i?|/o) and the bracket (/o|/o) are not defined. Thus, in principle, 
we cannot define the mean value of H on |/o) as we would do in ordinary Quantum Mechanics. This is precisely the 
sort of questions which we meant above. 

The aim of the present paper is to recall on the attempts to define the mean value of _ff on a Gamow state and 
discuss their advantages or disadvantages. These attempts use the following concepts: 

1. - The mean value of the energy of a decaying state must be zero, because the energy of a decaying process should 
be invariant. A non-zero energy will be in contradiction with the principle of conservation of the energy Jist p6| . 

2. - The energy average of a Gamow state should be complex because the Gamow state is an eigenvalue of H with 
complex energy. 

3. - If we admit that Gamow states are genuine quantum states they must have a real energy average that should 
be determined from first principles Q p7[ |. 

4. - Gamow states admit a representation as normalizable vectors in a Hilbert space in which the Hamiltonian is 
not a self adjoint operator and it has an spectrum of eigenvalues extending from —oo to oo. In this case, we can 
define scalar products of Gamow vectors and a mean value of the energy, which is real . 

In the last section, we shall present and evaluate these four possibilities. Before, in the next section, we recall the 
definition of a Gamow vector and enumerate some of its properties. 



II. THE GAMOW VECTORS AND THEIR PROPERTIES. 



In ordinary Quantum Mechanics, it is customary to associate a normalizable vector (with unitary norm ) to each 
pure quantum state. If we have a quantum decaying state and we want to assign to it a vector ?^ in a separable 
Hilbert space we can define its non-decay probability for t > as 

P{t) ^ where A{t) i:\e-'^"\i}; > . (1) 

Since ip is normalized one has < P{t) < 1 and the decay probability is given by 1 — P{t). A{t) is called the non-decay 
amplitude. It has been shown that A{t) is roughly exponential for values of time which are not too short neither too 
long. At t = the time derivative of A{t) is zero and this excludes the exponential time behaviour of ip for small 
times. For large values of time A{t) goes to zero slower than an exponential |p^ . 

This is what the theory predicts provided that a decaying state can be represented by a vector in a Hilbert space. We 
have already mentioned in the Introduction another argument to show that an exponentially decaying state cannot 
be represented by a vector in a Hilbert space. 

The experimental evidence on decaying states shows a decay which seems to be exponential up to the degree of 
experimental accuracy. In consequence, it is natural to consider exponentially decaying vector states as true physical 
states. But, we have to find out the mathematical nature of these objects. This problem has been satisfactorily solved 
by Arno Bohm and coworkers, in terms of rigged Hilbert spaces or Gelfand triplets |^ [?] ||2l| |^ [ p4| . 

In fact, we can look at a decaying process as the second half of a resonant scattering process in which the process of 
formation of a resonance is ignored Resonances can be defined in several ways, but it is generally accepted that, 
under rather general conditions, we can associate a resonance to each of the pairs of poles of the analytic continuation 
of the iS-matrix, S{E), on the real semi-axis of the energy values ^ [0 |]l8| |^^. In many implementable physical 
cases we can assume that in the resonant scattering process both the "free" Hamiltonian Hq and the "perturbed" 
Hamiltonian H = Hq -\- V, where ^ is a suitable potential, have the same continuous spectrum. This spectrum 
coincides with the positive semi-axis i?+ = [0, oo) and the identity of both continuous spectra is a consequence of the 
asymptotic completeness of the scattering. This property implies that the MoUer operators n± exist and fulfill the 
following relation |^] 

n±Ha = Hn± OT H ^ n±Hon^± (2) 

Furthermore, to simplify the formalism, we can assume that the energy spectrum is non-degenerate (which is the 
case of a spherically symmetric potential in the I = channel). In this case, it exists a unitary operator which 
connects the abstract Hilbert space of states and L^(i?+) (the Hilbert space of the energy representation) such that 
U HqU~^4){E) = E(f){E), where U diagonalizes the free Hamiltonian Hq. 

The explicit construction of the rigged Hilbert spaces for the decay process is relevant in our discussion and, 
therefore, we want to summarize it here. Further details can be found in the literature § (l| 111 H ||. First, 
we take the spaces of the so called very well behaved functions. These spaces are defined by analytic functions on the 
upper or the lower half planes of the complex plane C that vanish at infinity. The boundary values of these functions 
on the positive semi-axis i?^ are uniquely defined and, furthermore, it is possible to recover all the values of these 
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functions knowing their boundary values on i?+ ||2^. Let us call A± the spaces of these boundary values. We can 
construct a rigged Hilbert space suitable for the definition of decaying Gamow vector as follows. Firstly, we define 
$+ = fi+[/^^A+, where ft+ is the outgoing Moller operator. Then, if H is the Hilbert space of scattering states 0, 
we have the following rigged Hilbert space 

$+C?^C (3) 

Since every function in A+ determines uniquely an analytic function on the upper half plane for every vector 
(p+ E $+ we have an analytic function (j>{E)+ on the upper half plane. Its complex conjugate (j)'^{E) = is 
an analytic function on the lower half plane. Then, if = Eji — iT/2 denotes a resonance pole (F > 0) the Gamow 
vector l/o > can be defined as a functional on $+ a^ 

< (p+l/o >=0t(zfl) (4) 

This definition implies several important properties namely: (i) |/o >G is a continuous anti- linear functional 
on $+, (ii) l/o >^ 7i, (iii) the Hamiltonian H satisfies C $+ and is continuous with the topology on $+. It 

means that H can be continuously extended (with the weak topology) to a continuous operator on <I>^, so that the 
action of on |/o > is a well defined operation. Then, we have -ff|/o >= z^l/o >, which means that the Gamow 
vector l/o > is an eigenvector of H with eigenvalue z^. This is possible because <I>^ is not a Hilbert space, and, (iv) 
the time evolution of |/o > is only possible for positive values of t. We can show that M ||l3| 



e 



'**^|/o >=e-'*^«e-"|/o > (5) 



This means that |/o > decays exponentially. All these properties justify the choice of |/o > as a representation of the 
decaying Gamow vector. 

In a resonant scattering process, together with the decaying channel, it exists a process of capture or formation 
of a resonance. This is called the growing or capture process and it can be described by the evolution of another 
functional, which is the growing Gamow vector |/o >Q. This is a functional on $_ and therefore an element of the 
vector space As functional its definition coincides with the one given in Eq. (Q) by replacing the function (f)^{E) 
by (j)^{E), which is analytic on the lower half plane, and the point zji by its complex conjugate z'^ = Efj + iT/2 

<^Ah>^(^*{z*j,). (4') 

Its properties are the following: (i) |/o >^ 7i, (ii) H\Jq >— z'^jfo > and, (iii) its time evolution is well defined if i < 

i.e., l/o > increases exponentially until t = 0, which is conventionally the time at which the capture process is 
completed and the decay starts. 

From these definitions it can be concluded that Gamow vectors obey 

< fo\ip+ >^< ifi+lfo > and < 7o|(^_ (p_|/o >*, (7) 

where the star denotes complex conjugation. The bra vectors < /o| and < /o| are continuous linear functionals on 
$+ and respectively, 

< fo\H ^ z*j, < fo\; < fo\H ^ zn < fo\. (8) 
The time evolution of < /o| is defined for t > only and it gives 

</o|e"^=</o|e^*^-e-r* (9) 



■^I.e., the absolutely continuous space for the total Hamiltonian H. 
^Here, we are using the notation in 

*The question about the physical meaning of this functional i.e.: whether it is related to the capture or creation of a resonance, 
has not been answered yet. If capture and decaying processes are not equally probable they can not be symmetric in the sense 
presented here. Thus, in this interpretation, |/o > would be the time reversal of the Gamow vector |/o >. 
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and the time evolution of < /o| is defined for i < as 

<7o|e**^=</o|e'*^«er* (10) 

In this context, it seems that exponentially behaving state vectors are the only class of vectors which can represent 
Gamow states. To construct a representation, one has to consider these states and also a background^ physically 
produced by the interaction with the environment, re-scattering processes, etc p8| and mathematically by contour 
integrals in the complex plane |^0| |2|] |^^. Usually, the resonant scattering process is produced by the 

interaction of a free prepared state with a potential, creating the resonance. The prepared state must be represented 
by a Hilbert space vector, say ip. Then, if Q^ip = (/?_, we can show that if- is the sum of two contributions. One of 
these contributions is proportional to the decaying Gamow vector |/o > and the other to the background, represented 
by a vector in $+, so that flltl |l| |§ |l| ||] 

— a l/o > +|background > (11) 

where a is a complex number. This background would be responsible for the deviations of the exponential law on the 
range of short and large times. 

Our next goal is to present and compare the definitions of the energy average on Gamow vectors. 

III. DEFINITIONS OF ENERGY AVERAGE ON GAMOW VECTORS. 

In this section we shall review the known results obtained in dealing with the definition of energy averages on 
Gamow states, which are available in the literature, as well as our own definition of it. 

A. The mean value of the energy is equal to zero. 

It was Nakanishi who first proposed this idea [|l^. In fact, if H\fo >= zulfo > and < /o|-ff = z^j < /o|, we have 
that < fo\H\fo >= ZR < /ol/o >= 4, < /o|/o >• This implies that (zr - z}^) < /o|/o >= =^< /o|/o >= and, 
therefore, < /o|-ff|/o >= 0. The weak point of this argument is that the bracket < /o|/o > is not defined. Further 
attempts to define it have been made p8[ , but the results are not convincing from a mathematical point of view. 

B. The averages are complex. 

In specific models, like Friedrichs's model, the bracket < /o|/o > is well defined and its value is one ||2^. If we 
try to obtain this result in a general model independent setting we conclude that < /o|/o > can be defined as a 
distribution-kernel and that it has the value one, although it is not clear if this is the unique choice po| . If we now 
define 11 = |/o >< /o|, it is now obvious that 11^ = 11. This idempotency suggest us that that 11 could be taken as the 
density operator for the decaying Gamow vector |/o >• Now if it would be possible to define Tr {i?n} and this would 
be a candidate for the average value of on |/o >. In fact, with the help of some generalized spectral decompositions 
pO[ for the Hamiltonian in terms of the Gamow vectors and the generalized eigenvectors of H with eigenvalues in the 
continuous spectrum of H, we can define this trace in such a way that Tr{iJn} =< /o|-ff|/o > [0, thus 

< fo\H\fo >^ ZR < Mfo >= ZR (12) 

Yet this result cannot be acceptable from the physical point of view. Due to the time-energy uncertainty principle 
we cannot measure simultaneously the real part of zr, which is the resonant energy, and its imaginary part, which is 
proportional to the inverse of the half life. Thus, zr cannot be the average of any measurement process and cannot 
be accepted as the energy average. Also, from these considerations, we conclude that the energy average of a Gamow 
vector, if it can be defined, should be real. 
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C. The averages are real in the interpretation of Bohm. 



This point of view is based in the idea that it is possible to construct a rigged Hilbert space, in which the Gamow 
vector is a vector in the Hilbert space, under the following conditions: 

1. The continuous spectrum of H is the whole real axis, 

2. H is not self adjoint (although it is still symmetric, i.e., < (plHip >=< iJ^jV' > for all (f) and i/j in the domain of 
H), and, 

3. from the point of view of the Hilbert space, the Gamow vector is not in the domain of H, but the action of H 
on the Gamow vector is well defined in the dual space (that includes the Hilbert space). 

The spaces of analytic functions on a half plane that we are using here are spaces of Hardy functions ]3^ ] 
These functions are determined by their boundary values on the positive semi-axis — [0,oo). Moreover, we 
have an explicit formula to recover their values on the half plane (including the negative semi-axis, = (— oo,0]) 
from their values on the positive semi-axis If we denote by 7ij_, the spaces of Hardy functions on the upper (-I-) 
and lower (-) half planes, and by 7il|^+ the spaces of the restrictions of the functions of 7Yj_ on the positive semi-axis, 
there is a one to one mapping 9± ||7|] such that 

9±nl^nl\^^ (13) 

As a matter of fact, as our spaces of analytic functions, we take certain regular subspaces of H^. Let S be the space 
of all functions from R to C which are differentiable to all orders and that vanish at ±oo faster than the inverse of 
any polynomial (Schwartz space). Then, let us consider the spaces ^± = Hj. S. We have the following relation 
between and <f>± 

d± *± ' — > A±; V± = n±U-^ A± I — ^ $± , (14) 

or, equivalently, 

*± - V±e±^± (15) 

The spaces of Hardy functions Tij. are Hilbert spaces as subspaces of L^{R). Therefore, the norms in n\ and m 
L'^{R) coincide. The mappings 9± are one to one transformations from Hj. into L^(i?+). In this sense 9± are not 
unitary, from L'^{R) onto L'^{R^), because for any ip±{E) e Hj^ 



poo />oo 

\\0±v±\\mB+)= W±{E)\^dE< \^±[E)\^dE=\\^±\\L2(j,:^< 

Jo J ~oa 



(16) 



In fact, Lp± are boundary values of analytic functions and cannot be zero on the negative semi-axis unless they vanish 
identically. 

Now, we can construct a new rigged Hilbert space which is given by ^± C TYj. C The mappings 9±, induce 
two one to one mappings, 6*^, from onto A^, by means of the identity 

<0±v±\9lG±>=<^±\G±>, (17) 

where (p± G and G± £ 5*^. The mappings 6*^ are not extensions of 9± because of the non-unitary of 9±. On the 
other hand, the unitary of V± allows us to extend them into the dual spaces by means of a similar formula:, i.e. if 
4>± G A± then V±(j)± G ^± and if F± e A^ then V±F± G $± such that 

< V±^±\V±F± >=< 0±|F± > (18) 

Thus, to any (p G $± corresponds an analytic function 4>±{E) G and 4>±{E) = 9^^V^^ip±. Therefore, the 

Gamow vectors can be represented as vectors in ^± as {9^)^^V^^\fo > and {9^)^^VZ^\fo >■ However these 
formulas are unpractical to obtain the Gamow vectors in the new representation. In order to find them we shall use 
the definition of |/o > and |/o > and the Titchmarsh theorem |32j ||3^ |Q on Hardy functions. Then, we have 

1 r°° (h'^(F) 



and 
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< ^_|/o <l>t{zR) r $^dE (20) 



These formulas imply that 



and 



(''a-'v;-i/.>=5-^ (21) 

(,^,-.^-.|/„>._L_i_ (22) 

In this representation the Gamow vectors are square integrable, i.e., they belong to the Hilbert space L^{R). 
Therefore, we can define brackets and scalar products between them. We can also define energy averages on these 
vectors. This is, however, not an easy task. First of all, we must observe that in the new representation, the 
Hamiltonian H is given hy E — Vj^^HV±, where V± — V±9^^. It is easy to show that E(I}±{E) = E(j)±{E), i.e., the 

multiplication operator on L^{R). This is why we do not add subscripts in E. On Tij^, E is still symmetric but it is 
not self adjoint (has different deficiency indices on Tl^)- Its spectrum is purely continuous, simple and coincides with 
R. The definition of energy averages in this representation is hampered by the fact that 

E^— = — ^ (23) 

E-ZR E^ZR ' 

is not square integrable. This only means that the function {E — zr)~^ does not belong to the domain of the 
multiplication operator E. However, since E can be extended by continuity to the identity (|2^) makes sense. 

The representations (Eqs (^l|)) and (Eqs (22)) of the Gamow vectors, when properly normalized, could be used 
now to define a mean value of the energy for these states. The normalization that we are going to use is the Hilbert 
space normalization, i.e., if we call 

i^^^a-^; ^:<^^a-^, (24) 
h ~ zr £j - Zj^ 

then 

ll^'ll' = r 7^ - (25) 

J-oo {E - zr) + (r/2) 

Therefore, = 1 if a = 1/0?- 

Now, let us define the mean value of the energy on the decaying Gamow vector as 

<il]^\E\il^^ > (26) 

and let us evaluate its value. Since Ei/j-^{E) — E%Ij^{E), we have that 

< r\E\r >= - r — ^ dE^^ r (27) 

' ' TT J_^E~z*j^ E~ZR ttTJ' ^2 y ) 



\ r/2 ; 



The change of variables 



E — Er 

r/2 



transforms the last integral in (^7|) into 

Er /■°° dx r xdx 



(28) 



(29) 



TT 7-00 2^ +1 27r 

The first integral in (p9|) has the value it. The second admits a Cauchy principal value equal to zero. Thus, we find 
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<i^''\E\^j'' >=Er (30) 



and 

<i^'^\E\^'^ >^Er (31) 

We see that this definition of the energy average of Gamow vectors gives the same real value for both Gamow vectors 
and coincides with the resonant energy. In addition, due to the adopted normalization, we have that < i/;^!'!/;^ >=< 
^G|^G 1^ Furthermore 

<i:^\r>-- r -^ = (32) 
J-oo {E~ zrY 

Analogously, < >= 0. 



D. The averages are real in Berggren's interpretation. 

Berggren's approach to the mean value of the Hamiltonian on a Gamow state can be formulated very similarly 
to Bohm's. Following we shall not use Hardy functions to construct our Gelfand triplets. Instead, we consider 

here another triplet ^ C 7i C for which the space consists of tempered ultra-distributions. A simple definition 
of these objects can be found in Q |36j and a complete account in |^ [|8j. Vectors in ^ are entire analytic functions. 
Vectors in are represented by pairs of analytic functions on the open upper and lower half planes, respectively. If 
we call ipuiz) and 'ijji{z) these functions, we can write 

i:{E)^ME + iO)-ME-iO), (33) 

where ^u{E + iO) and 4'i{E — *0) represent the boundary limits of ^u{z) and ipiiz.) on the real axis, respectively. For 
any tp € S,^ and any z G C — i?, we have 

were we use in (^) the sign + or — for z on the upper or lower half plane, respectively. 

Now, we have a similar scheme to that presented in the previous section. The total Hamiltonian H is represented 
by the multiplication operator on ^ and it can be extended as a continuous operator into the dual The nuclear 
spectral theorem guarantees the existence of a complete set of eigenvectors \E > of H |^. As eigenvectors, they 
obey the identity H\E >= E\E >. The Gamow vectors are now defined a^ 

1 r \E > dE J 1 r \E > dE 

l/o >= -TT- / -E^ — — , l/o >= 77- / — -:r- (36) 



2771 7-00 E ^ ZR 2771 7-00 E — Zj^ 

Within this scheme, the mean value of the Gamow vectors is defined as in the previous section and it gives the same 
results]^. 

The coincidence between this last result for the mean value of the Gamow states and Bohm's one, presented in 
the last sub-section, comes from a re- interpretation of Berggren's definition given in In fact, for a spherically 

symmetric potential and for an arbitrary value of the angular momentum I, we can write the normalized decaying 
Gamow vector as 



^The decaying Gamow vector is denoted in [35-36] as \Eq > and the growing Gamow vector as \Eg >■ 

^In both cases, we can define the probability distribution associated to the Gamow states and it is given by (for the decaying 
Gamow vector) P{E) = | < E\fo > |^ = ^ ie-e^Y^+t'^ • same expression is found for the growing Gamow vector. 
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iA>=jEr>ir^*. (37, 

V TT Jo V m E[k) - ZR 

where, k = |k|, i?(k) = I? jlra and is the unit vector in the direction of k. For the growing Gamow vector, we have 

~ f2r IT \k,k,i> „ 

/o>--m/— / J- ' , dk. 38 

Now, let A be an arbitrary observable. We can define the mean value of A on |/o > as 

^i.i^ p^, /•°°^,,VfcF <k',k',i'\A\k,k,i> 

^ Jjf Ja Jo ^ {E{\<.') - ZR){E{]<i) - z*j^) 

If we replace A by 7J we obtain, straightforwardly, the value Er for this average. Instead, Berggren defines 
this mean value as Real < /o|A|/o >. As a matter of fact, one can easily show that 

< h\A\h >- Real{< /o|A|/o >} + o{T^), (40) 
which means that Berggren's approximation coincides with Bohm's to the first order in F. 



IV. CONCLUSIONS 



In this paper we have compared different definitions of Gamow vectors. We have shown the equivalence between 
Bohm's and Berggren's definitions of the mean value of the Hamiltonian on a resonant state. Our main result, 
concerning this equivalence, is the realization of the average value of the Hamiltonian on a resonance as a real 
function which depends on both the real and the imaginary parts of the complex energy. This result is supported, 
mathematically, by a proper treatment of Gamow vectors in a rigged Hilbert space. 
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